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Abstract We observe that the Laplacian of a random graph Gn on N vertices represents 
an explicitly solvable model in the limit of infinitely increasing N . Namely, we derive 
recurrent relations for the limiting averaged moments of the adjacency matrix of Gm when 
N — > oo. These relations allow one to study the corresponding eigenvalue distribution 
function; we show that its density has a infinite support in contrast to the case of the 
£>. . ordinary discrete Laplacian. 
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The spectral theory of graphs attracts more and more attention in mathe- 
matical physics (see e.g. |], 0). It is widely believed that the spectral prop- 
Q_i' erties of graphs are related with their geometry that is an intriguing assertion 

from various points of view. However in general, without the knowledge of the 
graph structure, one can prove just abstract theorems with not large number 
of concrete results. In this connection, the study of the spectrum of random 
graphs, as the typical in certain sense objects, can be useful. We confirm this 
assumption by showing that the models related with large random graphs are 
rS \ explicitly solvable. 

Given a set Vn of N vertices v%, . . . Vn, one can assume that some pairs of 
them are either joined by one undirected edge or are not joined. Then Vn and 
the set E N of existing edges determines the graph r^ = (Vat, E n ). In this case 
the graph Laplacian A^ acting on functions f(v),v G Vn is determined by the 
formula 

A ff /K)= E [/(**)-/(«*)], (!) 

where the sum runs over the vertices joined with V{. 

It is clear that An has N x N real symmetric matrix and 
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y(N) _ A (N) 



where A^ is the adjacency matrix of the graph T 
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%1 \ 0, otherwise, 

and V^ is the diagonal matrix such that V- = deg(fj) is the number of 
edges attached to t>j. Note that A u = and that 

A N = P N [A^} 2 P N -A^\ (3) 

where P/v is the orthogonal projection onto the diagonal. Often the set of 
eigenvalues Ai < . . . < Ajv of A N is called the graph's spectrum ||. Firstly 
we concentrate on this problem and turn to the spectrum of A^ at the end of 
this article. 

Our principal goal is to describe results obtained for the adjacency matrix 
A( n >p) of a random graph G% = (Vjv, &$ (u;)). We are related with the well- 
known model (see e.g. Q), where {G% } represents the ensemble of graphs 
where each edge, independently from the others, is either present or absent 
with probabilities p/N, < p < N and 1 — p/N, respectively. Thus we get a 
real symmetric matrix A^ N ' P ^ with the entries {A\, , i < j} determined as a 
family of jointly independent random variables such that 

A<> N <ri = n nT1 H A^ N ^ - / 1 » with Probability p/N, ,.. 

u v ' 1 0, with probability 1 - p/N. y ' 

We are interested in asymptotic behaviour in the limit N — > oo of the moments 
M (N,p) = Ejjy-i Tr A ay } A = A {N, P )^ ^ 

where 

1 N 
N -i Tr [A (N, P)]S = £ AxiX2 A X2X3 ■ ■ ■ A Xa _ lX1 . (5b) 

Xi = l 

We consider for simplicity the case of p = 1. Our main result is that 

n^oo s 10, ifs = 2A; + l 

and rrik can be found from equality 

k 

m k = J2 W k(r) (6a) 

r=0 



where the numbers W&(r), k > r > are given by the following recurrent 
relations 

^) = EE E w«-*-i(«) )_ '^-0, (66) 

i=\ j=v—i 1=0 \ / \ / 

supplied with the initial condition Wj(0) = 5j . 

Derivation of (6) essentially uses the fact that A^ N ' P "> (4) is the matrix whose 
entries are jointly independent (excepting the symmetry condition) random 
variables. For such class of matrices, the limits of moments (5) were studied 
first by E. Wigner ||. The ensemble considered in || is somewhat different 
from (4), and the method proposed there and then developed in is applicable 
to (4) in the asymptotic regime p ^> 1 only. 

The principal observation of || and [[J is that the sum over the sequences 
X 2 k = (xi, X2, ■ ■ ■ , x 2 k~i, X\) in (5a) can be replaced by the sum over the equiv- 
alence classes of the walks X 2 k', these classes are labelled by the plane rooted 
trees r G T k drawn in the upper half-plane with the help of k edges. 

In paper the method of || and was further improved up to the 
stage when p ~ 1 can be also included into consideration. Summimg up the 
arguments of P, || and 0, we claim that 

fc— l 

™* = E E u(k]T) f (?) 

q=0 rdT k _ q 

where U(k; r) is the number of all possible different walks of 2k steps such that 
cover the tree r G T k - q according to the following 

Riding rule: the walk starts and ends at the root and each edge of the tree 
is passed even positive number of times; when there is a choice where to move 
further, the walk chooses either one of the edges that has been already passed 
or the most left edge among those that have not yet been passed. 

Let us call the walks of this type the even ordered walks covering the tree. 

Demonstration of (6). 

It is convenient to imagine that some particle moves in the upper half-plane. 
It starts at the root p, then makes 2k steps and returns to the root; each step 
either produces a new vertex or ends at the vertex that already exists and is 
connected by the present vertex by an edge. The condition is that each step 
made in one direction is repeated in inverse direction. 

It is easy to see that there exists a tree r G T k _ q (and the only one) that the 



walk of the particle represents one of the coverings of r satisfying the riding 
rule. 

Let us introduce the set W u (v ) of all possible walks of 2u steps that return 
to the root v times and denote by W u (v) = \W u (v)\ the number of such walks. 
With this definition in mind, it is not hard to derive (6) from (7). 

Indeed, regarding the set of walks W u (v ), let us denote their first step by 
(p, a) and consider those walks that pass this edge 2i > 2 times. It is easy 
to see that a by itself can be regarded as the root of a subwalk that starts 
and ends at a and never visits p. Chronologically, this trajectory consists of 
possibly several pieces. Let us call their union the first subwalk. Similarly, the 
second subwalk starts and ends at p and never visits a. 

Let us assume that the second subwalk belongs to Wj(v — i) and the first 
subwalk belongs to W u -i-j(l)- Then (6b) follows, where in the product, the 
second factor represents the number of possibilities to perform i — 1 steps from 
a to p under the condition that we arrive I times to a by the first subwalk; 
the last passage a — > p is made after that the first subwalk is completed. The 
third factor in (6b) corresponds to i — 1 steps p —>■ a made after v — i arrivals 
to p by the second subwalk. □ 

In the following remarks we briefly discuss the result obtained 

1. The first conclusion relates the maximal eigenvalue of A^' 1 ' . It is clear 
that the number of non-zero elements in a row is asymptotically Poissonian 
when N — > oo. Thus, with non-zero probability, in (Viv,i% ) there exists 
a vertex that has large number of edges attached to it. This explains why 
H^W 1 )]! — > oo (see also 0) while the spectrum of the ordinary discrete Lapla- 
cians (or more generally, of the regular graphs) remains bounded. 

The stress of the present article is that there is an infinite number of the 
eigenvalues of A( N,P ' that tend to infinity as iV — » oo. We derive this conclusion 
from relations (6) and the representation 

M Wp) = E | AMajVp(A)) (8) 

where aiy p (\) = cr(X,A^ N ' p ') is the normalized eigenvalue counting function 
«r(A) = #{A; < \}N~\ 

Basing on (6b), one can easily show that Wk(r) < (C\k) 2r where C\ is some 
constant. Then we derive from (6a) the estimate 

m 2k < (C 2 k) 2k 

that implies the uniqueness of the limiting measure a^ = lim^r cr^ 1 . 



^From another hand, it is not hard to derive the estimate 

m 2k > (k/2) k/2 (l + o(l)), k -, oo. (9) 

Indeed, let us take the set Wk(k) and consider the walks that cover one par- 
ticular tree r' k , 2 of k/2 edges, namely that has all k/2 edges attached to the 
root (here we assume k — 21 for simplicity). Among these walks there are 
those that pass each edge four times. It clear that the number of these walks 
is 0/2)! that results in (9). 

Inequality (9) implies that the measure da'- 1 '' has an infinite support. This 
means that the mathematical expectation of the number of eigenvalues that 

go to infinity is of the order 0(N) as iV — > oo. 

(p) 

2. Let us note that one can easily apply our method to study the limits rn£ 

of moments (8) with arbitrary constant number p. In paper || the arguments 
have been presented showing that there exists a critical value p c ~ 2 that sepa- 
rates two different pictures of the limiting spectra of A^ N ' P \ This is supported 
by the facts from the random graphs theory stating that p' c = 1 is critical for 
existing of an infinite connected component in (Vat, En ) as N — > oo (see e.g. 
@]). Our preliminary numerical studies of the moments m k show that a^ is 
not sensitive up to the level that reflects the role of p c or p' c . 

3. It should be mentioned that the moments (7) we compute represent a 
special kind of the matrix integrals. If one considers the matrices A^ (4) as 
having the Gaussian invariant distribution, then the corresponding moments 
(5) can be computed explicitly for fixed N and the formulas obtained are 
related with the number of maps drawn on the surfaces of given genus (see e.g. 
for an accessible introduction). 

In this aspect, our model characterises another class of random matrices 
with independent entries that have explicitly countable integrals over them 
(see |TI| for related recent results concerning their asymptotic properties). 

4. Finally, let us indicate the way to find explicit relations for the moments 
of the random graph's Laplacian. Regarding definition (3), we see that one 
faces the problem of computing the mixed moments 

E {A r ^ (xi, x 2 ) [A 2 (x 2 , x 2 )] Sl • ■ • A r > ( Xj , x x ) [A 2 (x h xt)} 8 * } . 

One can apply again the method of |5], |6| with the modification that the trees 
are constructed with the help of edges of two different colours and the edges of 
one of the class can be used as the leaves only. Certainly, in the case of finite 



p the computations become more difficult than those of (6), but the general 
approach of ]7J and the present article remains valid. 

Note that the case of multiple edges can be also included into consideration 
by allowing Aij (4) to take positive integer values. Moreover, our method still 
works for the random graphs with independent random weigths on edges; then 
the adjacency matrix (2) takes the form B^ = EyAy, where Ey,i < j are 
jointly independent random variables. All these questions will be considered 
in more details in separate publication. 
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